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Abstract 

Within the framework of generalized collinear factorization in perturbative QCD (pQCD), we 
study the effect of initial multiple parton scattering and induced parton energy loss in Drell- 
Yan (DY) process in proton- nucleus collisions. We express the contribution from multiple parton 
scattering and induced gluon radiation to the DY dilepton spectra in terms of nuclear modified 
effective beam quark distribution functions. The modification depends on the quark transport 
parameter in nuclear medium. This is similar to the final-state multiple parton scattering in 
deeply inelastic scattering (DIS) off large nuclei and leads to the suppression of the Drell-Yan 
cross section m. p -\- A relative io p + p collisions. With the value of quark transport parameter 
determined from the nuclear modification of single-inclusive DIS hadron spectra as measured by the 
HERMES experiment, we calculate DY spectra \n.p + A collisions and find the nuclear suppression 
due to beam parton energy loss negligible at the Fermilab energy Siab=800 GeV in the kinematic 
region as covered by the E866 experiment. Most of the observed nuclear suppression of DY spectra 
in E866 experiment can be described well by parton shadowing in target nuclei as given by the 
EPS08 parameterization. The effect of beam parton energy loss, however, becomes significant for 
DY lepton pairs with large beam parton momentum fraction x' or small target parton momentum 
fraction x. We also predict the DY cross section m. p + A collisions at lower beam proton energy 
-Eiab=120 GeV and show significant suppression due to initial state parton energy loss at moderately 
large x' where the effect of parton shadowing is very small. 
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I. INTRODUCTION 



A basic assumption in the coUinear factorized parton model in perturbative QCD (pQCD) 
is the small intrinsic transverse momentum of initial partons inside beam hadrons or nuclei 
relative to both partons' longitudinal momenta and the large energy-momentum scale Q 
involved in hard partonic interactions. To go beyond such a coUinear factorized pQCD model, 
one has to expand the hard partonic part of the hard interaction in the intrinsic parton 
transverse momentum. In such an expansion, the first term gives rise to the normal coUinear 
factorized pQCD results, known as the leading twist contributions which will depend on 
leading twist transverse-momentum-integrated parton distribution functions. The higher 
order terms in the Taylor expansion can be combined together with contributions from 
multiple parton interaction between incoming or out-going partons and the remanent of 
target hadrons/nuclei. These contributions are known as higher twist contributions which 
can be expressed as the convolution of hard parts involving multiple parton scattering and 
multiple parton correlation functions inside the target hadrons/nuclei in the framework 
known as generalized collinear parton model [H, [2|. These higher- twist contributions in 
general are suppressed by powers of the momentum scale in the hard processes, 

In hard processes involving a large nuclear target, higher-twist contributions often depend 
on multiple parton correlation functions inside the target nucleus which are enhanced by 
the nuclear size Ra ~ A^/^ These will give rise to a leading and non-trivial nuclear 

dependance of the higher-twist contributions to the hard processes beyond the incoherent 
superposition of nucleons inside a nucleus. One example of such nuclear enhanced higher- 
twist contributions is the suppression of leading hadron spectra in deeply inelastic scattering 
(DIS) off a nucleus relative to that in DIS off a nucleon and the suppression is observed 
to increase with the nuclear size as in the HERMES experiment JgII. Such a suppression 
has been studied within the generalized collinear parton model [7H9| and was shown to 
be caused by multiple parton scattering of the struck quark and induced gluon radiation 
inside the target nucleus. The gluon radiation induced by multiple parton scattering that 
reduces the energy of the propagating quark can be expressed as a nuclear modification 
to the effective fragmentation function of the struck quark [7], which is determined by the 
quark-gluon correlation function inside the nucleus, a higher twist matrix element of the 
nucleus not calculable in pQCD. Assuming a large and loosely bound nucleus, one can 
approximate these higher twist matrix elements as a product of initial quark distribution 
and gluon density distribution inside the nucleus, the latter is also known as the quark 
transport parameter inside the nuclear medium The parton transport parameter also 
describes the medium broadening of transverse momentum squared [lO, 11| per unit distance 



and is equivalent to the saturation scale per unit length in a large nucleus [12|, Il3[. One 
can determine this quark transport parameter by comparing the leading hadron spectra 
suppression in DIS off a nuclear target to experimental data. The predicted nuclear size and 
quark energy dependence of the hadron suppression agree very well with the experiment 
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One can further extend the study of medium modification of parton fragmentation func- 
tions and parton energy loss or jet quenching due to multiple parton scattering and induced 
gluon radiation to the case of parton propagation inside a hot quark-gluon plasma (QGP) 
in heavy-ion collisions. In this case one has to replace the cold nuclei with a hot thermal 
medium of quarks and gluons j^, 14, 16]. The parton transport parameter is then propor- 
tional to gluon distribution density inside the hot QGP. Therefore the study of jet quenching 
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in high-energy heavy-ion colhsions can provide crucial information of s pac e-time profile of 



parton density distribution in the early stage of heavy-ion collisions [17(|. The observed 
suppression of single [isl, 19 1, dihadron spectra j20| and gamma-hadron correlation j2l|, 22| 
in high-energy heavy- ion collisions at RHIC are all consequences of strong jet quenching 
and provided the evidence for the formation of a strongly coupled QGP in the center of 
heavy- ion collisions at RHIC 2^ 24 1. Recent data from heavy-ion collisions at the LHC [Isf 
have also shown a wide range of phenomena of jet quenching which point to the formation 
of a strongly coupled QGP similar to that formed in heavy-ion collisions at RHIC. 

Multiple scattering and induced parton energy loss have been studied in several theo- 
retical frameworks [26l^29|. among which the higher-twist approach can be applied to jet 



propagation in both hot QGP and cold nuclei. In high-energy heavy-ion collisions, ener- 
getic partons which one uses as hard probes have to traverse not only the hot QGP in the 
early stage of the dense matter but also hadronic medium after hadronization of the QGP 
matter. One therefore has to include the effect of jet quenching due to multiple parton 
scattering during the hadronic phase of the dense matter evolution in high-energy heavy- 
ion collisions. Using the information of parton energy loss in cold nuclei in DIS off large 
nuclei and extrapolate to hadronic medium at finite temperature, a recent study [soj finds 
non- negligible contribution from jet quenching in hadronic phase to the final suppression of 
large pt hadron spectra in heavy-ion collisions at RHIC and LHC. Therefore, the study and 
inclusion of parton energy loss in cold nuclei and hadronic matter should be an integral part 
of phenomenological study of jet quenching in heavy-ion collisions. 

In this paper, we will study the effect of multiple parton scattering and parton energy loss 
in the Drell-Yan (DY) process of dilepton production in p + A collisions. Since dilepton pairs 
in the final state do not have strong interaction with the target nucleus, only initial state 
interactions are important in DY process. Therefore, it is an ideal tool to study fast parton 
propagation and energy loss in cold nuclear medium 3lH37| . The higher-twist approach has 
also been applied to the study of nuclear effects in the DY process in p + A collisions jsS, 39 . 
These studies were restricted to the large transverse momentum DY pairs, where one can 
neglect the interferences associated with induced gluon radiation. In the calculation of the 
total DY cross section, one has to integrate over the transverse momentum of the DY pairs 
and therefore has to take into account of the interferences. 

Our approach to the nuclear effects in DY pair production in p+A collisions is very similar 
to the study of nuclear modification of the effective fragmentation functions in DIS. Initial 
multiple parton scattering and induced gluon radiation before the DY pair production via 
quark-anti-quark annihilation can lead to parton energy loss and attenuation of the incoming 
beam partons. In the collinear approximation, transverse momentum of the initial partons 
are neglected. The transverse momentum of the final DY pairs is then determined by the 
transverse momentum of the radiated gluon either in the vacuum or medium induced gluon 
bremsstrahlung. In the process involving multiple parton scattering, there are destructive 
interferences between the amplitudes of vacuum gluon bremsstrahlung and gluon radiation 
induced by multiple scattering. It is often referred to as the Landau-Pomercanchuk-Migdal 
(LPM) interference |4o[ . The LPM interference effect can be characterized by the formation 
time of the radiated gluon which is inversely proportional to the transverse momentum 
squared For large or short formation time, interference between vacuum and medium 
induced gluon bremsstrahlung is negligible. However, for small q^, the gluon formation 
time becomes longer than or comparable to the nuclear size and there is strong destructive 
interference between vacuum and medium induced gluon bremsstrahlung that cancels the 



3 



effect of gluon radiation due to multiple parton scattering. In our study of DY pair spectra, 
we will include explicitly the LPM interference within the framework higher-twist approach. 

Similar to the parton energy loss and medium modification of the effective parton frag- 
mentation functions in DIS, we can express the nuclear effect in the DY process due to the 
initial multiple parton scattering and induced parton energy loss by the beam parton in 
terms of modified beam parton distribution functions. The medium modification can be 
effectively included in terms of the modified splitting functions in the QCD evolution equa- 
tions which depend on twist-4 quark-gluon correlation functions inside the nucleus. With 
the approximation that the quark-gluon parton correlation can be factorized as a product of 
the quark and gluon distribution, the medium modified splitting functions depend only on 
the parton transport parameter g, which is related to the distribution density of soft gluons 
inside the nucleus or the broadening of transverse momentum squared per unit length due 
to multiple scattering. This is the only free parameter in the higher-twist study in both 
DIS and DY process with a nucleus target. One therefore can use the information from 



the phenomenological study of hadron suppression in DIS experiment [15| and calculate the 



nuclear suppression due to parton energy loss in the DY process [41 . 

Coherence in multiple scattering between beam and target partons can also lead to the 
nuclear shadowing of parton distributions inside the target nucleus. This and other nu- 
clear modifications of the parton distributions has been measured in nuclear DIS and 



parametrized [43|, |4^ from global fitting of experimental data involving nuclear targets. 



These nuclear modifications of parton distributions can also lead to nuclear effects in DY 



spectra in p + A versus that in p + p collisions. We will use the EPS08 parameterization [43 
of the nuclear parton distribution functions (nPDF) for the target nucleus in the calculation 
of DY spectra in addition to the effect of parton energy loss and medium modification of 
the beam parton distributions. We will then compare the calculated DY spectra to the data 
from the Fermilab Experiment 866 (E866) [32]. We will compare the effect of nuclear shad- 
owing of target partons and the beam parton energy loss and investigate in which kinetic 
region the effect of beam parton energy loss dominate. 

The remainder of the paper is organized as the following. In Sec. [Tll we define our no- 
tations and discuss the leading twist processes and the vacuum Dokshitzer-Gribov-Lipatov- 



Altarelli-Parisi (DGLAP) evolution 45|, |46| of parton distribution functions. In Sec. IIIH we 
calculate the twist-4 contribution to DY cross section within the framework of generalized 
factorization. We then define the effective modified beam quark distribution functions in 
Sec. IIVI and discuss the energy loss effect on the DY cross section. We summarize our 
results in Sec. |Vl The complete results of the twist-4 annihilation-like and Compton-like 
processes are given in the Appendix |X] and [HI respectively. We also demonstrate that the 
twist-4 contributions from double-hard scattering can be understood as the successive single 
scatterings in Appendix [O 

II. LEADING TWIST CONTRIBUTIONS 

In this paper, we focus on the production of DY lepton pairs in unpolarized h + A 
collisions, 

h{p') + A{p) ^ ^' (1) 
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whose cross section can be expressed in terms of the virtual photon 7*(g) production cross 
section, 

Here p' is the four-momentum of the incoming beam hadron, p is the momentum per nucleon 
inside the target nucleus and q is the four-momentum of the DY lepton pair. In the center- 
of-mass frame, we can choose the nucleus and the beam hadron moving along the "plus" and 
"minus" direction, respectively. Using light-cone notations, the momenta of beam hadron, 
target nucleus and virtual photon can be expressed as 

p' = [o,p'-,oj, p= [p+,o,oj, q=[{Q^ + qT)/'^q',q',qT], (3) 

where the "plus" and "minus" components are defined as p^ = (po ^Pz)/V^- 



According to the coUinear factorization |47|, the Drell-Yan cross section in h + A collisions 
can be expressed, 

%f -i:/<^^7.,.M%f ./<;.'^^. (4) 

in terms of the product of parton(c)-nucleus cross section dacA^-y and the normal twist-two 
beam parton distribution function fc/h{x'), where x' = q~ /p'^ is the momentum fraction 
carried by the annihilated beam quark. 

Inside the nuclear target, the beam parton may undergo additional scattering before the 
production of the lepton pair. Within the framework of twist expansion, one can calculate 
the contribution to the parton-nucleus cross section as a sum of contributions from single 
scattering, double scattering and even higher multiple scattering 

d^q ~ d^q d^q ^ ^ 

We will calculate contributions up to double scatterings in this paper. 




FIG. 1: Lowest order Drell-Yan process at leading twist. 



The leading-twist contribution to DY dilepton pair production at the lowest order C(a;°) 
comes from the quark anti-quark annihilation process as illustrated in Fig. [1] The corre- 
sponding differential cross section is given by |48| 



dQ'^dx' 



Yl j dxfg/A{x)fg/h{x')Ho{x,p,q) , (6) 
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where the summation is over all possible quark and anti-quark and fq/A^x) denotes the anti- 
quark (quark) distribution function in the nucleus with momentum fraction x and the hard 
part of quark anti-quark annihilation is 

TT / \ 2'/TG;em6„ 

Ho[x,p,q) = ^— 6{x-Xb). (7) 

6X'S 

Here xb = I'^V ' Q is the Bjorken variable in the DY process, s = {p + p'Y is the center- 
of-mass energy per h + N collision. 

To consider the scale dependence of beam quark distribution function, one has to con- 
sider radiative corrections to the annihilation processes and Compton scattering process as 
illustrated in Fig. In the annihilation process, the beam quark can radiate a gluon before 
its annihilation with the anti-quark from the target as illustrated in Fig. M^a) and (b). The 
annihilating beam quark can also be generated from a beam gluon splitting as shown in Fig. 
[2t^c) which is actually a Compton scattering process. 




{a) (b) (c) 



FIG. 2: Next-to- leading order Drell-Yan annihilation process (a,b) and compton scattering 
(c) at leading twist. 



After including the radiative contributions in Fig. [2] with the leading logarithmic approx- 
imation (for up to a factorization scale /x^), the renormalized beam quark distribution 
can be defined as 



fq/h{x',fi ) = fq/h{x') + 



as f^'^ dql 



271 Jo qq- Jx' C, 



fq/hiOlq^qgi^'/O 



(8) 



The first and second term inside the square brackets represent contributions from the anni- 
hilation and Compton processes, respectively. The corresponding quark and gluon splitting 
functions are given by 



(9) 



and 



^yg^qqiz) = Tr + {1 - z) 



(10) 



where Ci? = 4/3 and Tr = 1/2 are the color factors. 
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The renormalized quark distribution function satisfies tlie vacuum DGLAP evolution 
equation, 



ain/x2 2ttJx' ^ 



(11) 



Similarly, one can also obtain the renormalized target quark distribution function and its 
vacuum DGLAP evolution equations. 



III. TWIST-FOUR CONTRIBUTION 

When a parton is passing through nuclear matter, it will suffer multiple scattering and 
parton energy loss. In principle, contributions from multiple scattering are power-suppressed 
for large invariant mass M = Q or transverse momentum qt of the DY lepton pair. However, 
if the involved target partons in multiple scattering come from different nucleons inside the 
nucleus, some contributions can be enhanced by the nuclear size A^/^ for each additional 
scattering. In this paper, we consider twist-4 contributions which are enhanced by the 
nuclear medium size. We will only consider the case of secondary scattering with another 
target gluon because the involved gluon-quark correlation in the nucleus is believed to be 



much larger than quark-quark correlation in the double quark rescattering [49 . 



A. Factorization 




In the covariant gauge, the double scattering as shown in Fig. [3] in general can be cast 
in the following factorized form: 

d^q ~ 2eJ 2vr 27r 27r (27r)2^ , Z/i , Z/2 , fcT,P, gj 

X l(A|A+(|/2-,0T)V^,(0)7+V^,(r)^+(z/r,Z/T)|A), (12) 

where ko = [{x + X2)p^ ,0,0±], k = [xp~^ ,0,0±\, and we keep for now the relative intrinsic 
transverse momentum kx carried by target gluon with ki = [xip'^, 0, kr] and k2 = [{xi — 
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X2)p~^, 0, kj]- H is the Fourier transform of the hard part H in momentum space, 



H{y ,yi ,y2,kT,p,q) 



dxdxidx2e'''P y e'^^'P 



X H{^,x,Xi,X2,kT,p,p,q). 



-i{xi-X2)p+y2 



(13) 

As part of the twist expansion, one makes a Taylor expansion of the partonic part H in 
the gluon's intrinsic transverse momentum around kx = j50| : 



Hiy ,y^,y2,kT,p,q) = H{y ,y^ ,y2 ,kT = 0,p,q) 



+ 



OH 



dk^ 



krp -\- — 



dm 



2 dk^dkt, 



f\jrpf\jrj-} ~|~ ... . 



(14) 



fcT=0 



kT=0 urvrpiyrvrp 

The first term on the right-hand side of the above equation can be included in the leading- 
twist eikonal contribution as the gauge link, which makes twist-two parton matrix element 
from the single-scattering gauge invariant. The second term vanishes after integrating over 
kT for unpolarized initial and final states. The third term will give a finite contribution to 
the double-scattering, which, after integrating over d'^kx, becomes 

dy- dyi dy2 1 



d'^q 



1 



X [--9 



Q/3 



27r 27r 2ti 2 

"1 



2 dki^dk'. 



(A|F+(y2-)^,(0)7+^,(2/-)F+"(2/r)l^) 
'H{y',yi,y2,kT,p,q) 



(15) 



where the factor k^k!^ has been converted into spatial derivatives in the the field strength 
tensor after partial integration. 



B. Double scattering: Annihilation- like and compton-like processes 

There are two distinctive processes at twist-four which we will consider here. One is the 
annihilation-like processes as shown in Fig. HJ Before annihilation into a virtual photon, 
the beam quark may scatter with another gluon from the nucleus. Such multiple parton 
scattering processes will effectively modify the distribution of the beam quark. In our 
calculation, we neglect the leading order twist-4 processes because our purpose in this paper 
is to study the effect of radiative parton energy loss due to the multiple scattering on 
the beam quark distribution. Without induced gluon radiation, the leading order double 
scattering mainly contributes to the transverse momentum broadening of the leading parton 

As an illustration of the calculation involved, we consider the contribution from the first 
diagram in Fig. |H Using the Feynman rules, one can write down the hard part of this 
diagram 

Haii = J cia;da;irfa;2e'"P''^"e*"^P''^i"e-'("^-"^)^'^2/2-lT^[^^.y^^^^]r(2), (16) 



where. 



TT Cf 2 2 7 • {jp' + k2) g l-jq-ko) ^1 

2N^ {^P + ^2) — 1£ [q — kq) — 2 
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FIG. 4: Next-to- leading order Drell-Yan annihilation process at twist-four. 



Throughout this paper, we use the subscript "Aij" to denote different central-cut diagrams 
at twist-4. We have the following conventions: "A" stands for the annihilation-like processes 
while "i" and "j" correspond to the i-th and j-th diagrams shown in Fig. H] which are the left 
and right side of the cut diagrams, respectively. Similarly, for the Compton-like processes, 
we use the subscript "Cij" to distinguish different central-cut diagrams. 
For central-cut diagrams, the final state phase space is 



r(2) 



Xb 



(27r)3Q2 I _ ^ 
The momentum fractions are defined as 



5{x + Xi — Xb — Xt — xc)- 



(18) 



Xb Qt 



Xb z{k^ — 2qT ■ kx) 



xd 



Xb ,2 



Performing the contour integration / dxdxidx2, the hard part becomes 



(19) 



Hau 

I All 



c 



1 



-I All 



e{-y2)0{y- 



Vi 



-\pi{xB+xt+xc~XD)v'^y ixDV'^iy-i -Vi ) 



X 



-i{xt+xc-xn)p^{y -Vi ) 



-i(xt+xc-XD)p^y2 



(20) 



(21) 



There are four terms in the above result which correspond to the contributions from four 
possible physical processes, i.e., soft-hard, double-hard and their interferences. This comes 
out naturally because of the four possible choices of poles when performing the contour in- 
tegrations. The interference are destructive in small qt region, which cancels contributions 
from soft-hard and double-hard processes. In the limit — )■ 0, there is a complete cance- 
lation among all the processes. This destructive effect is important in our calculation since 
we are interested in the transverse momentum integrated spectra to which the small angle 
gluon radiation also contributes. 

According to the scheme in the twist expansion, the twist-four contributions comes from 
the second derivatives of H with respect to kx or the quadratic terms in the Taylor expansion. 
In general, one can rewrite the quadratic term in collinear expansion as 



dk^dk^n 



dxl^ dk^ dk^ 



H + 



dl d'^Xkj, 
dxkT dk^dk^ 



H + 



dl dx 



dXkrj 



dk^ 



OH 



+ 



dl dxkj 
dxkT dk^ J 



dH 



(22) 



9 



where H is the truncated hard part of the scattering matrix element when / that contains 
all the phase factors is factored out. We denote x^j, the dependent fractional momentum, 
such as xciXd in Eq. (!2T|) . As we will show later, the last term of the above equation that 
contains second derivative of the truncated hard part H with respect to the intrinsic trans- 
verse momentum is proportional to a quark-gluon correlation matrix element. Other 
terms in the above equations that contain derivatives of the phase factor / will be propor- 
tional to the derivative of the quark-gluon matrix elements with respective to the partons' 
momentum fraction. These terms are also proportional to the power of the momentum frac- 
tion xt ~ (ItIQ"^ &^c1 therefore are suppressed relative to the last term in the above equation 
for small Qt/Q"^ ^ 1- This feature is well illustrated by the the second-order derivative of 
Hah, 



-4» 



H 



All 



kT=0 



^2 2^F i 



1 



z 



I ^^ \ ^^All , 2 2 

-|- z[l — z)xt— V z x^ 



z^Iau + z(l - 2z)xt 



dl 



All 



dxc 



dx 



D 



d^Aii 

dxl; 



(23) 



xc=0,xo=0 



of which the leading contribution in the limit a;^ ~ ^ 1 comes from the first term. 

However, in the limit Xtp'^ <^ ^/Rai LPM interference becomes important and all the terms 
are equally important. In this region, non-perturbative effects also become dominant and the 
medium modification cannot be calculable anymore. In our higher-twist approach, we simply 
assume the medium modification is still given by the first term in the above equation that 
do not involve derivatives of the quark-gluon correlation functions. Under this assumption, 
the hard part from the first diagram in Fig. H] is then 



H 



All 



1 c 

4^ dk^dk^ 



Up 



z 



e{-y2)e{y- 



Vi je 



-ixtp+{y -y^ ) 



(24) 

Calculation for the rest of the diagrams and their interferences in Fig. H] can be similarly 
carried out. Summing these contributions together, one can obtain the twist-4 contributions 
from the annihilation-like processes. 



1 

-^9 



4"" dk^dk^T 



«e.e^«^i±^^^(-y,-)^(r - yr)e^^-+-)-"^" 



kT=0 



l-z 



Cpz^ (l - e-*^*^'''^^ )) (l - e~"'P'^y^ ) 



_j_ (jpQ-ixtp+{y -y^ )^-ixtp+y2 

Ca- 
2 



+ 0{ql/Q') 



-ixtp^y^ 



(25) 



where the color factor Ca = 3. The first two terms correspond to the first two diagrams in 
Fig. m respectively. The interferences between these two diagrams are given by the third 
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and fourth term. However, the diagram containing triple-gluon interaction gives power 
suppressed contributions Qt/Q"^ as compared to the other two diagrams and is therefore 
neglected. This is quite different from the final state multiple interaction in the nuclear 
DIS |3] where the dominant contribution due to induced gluon radiation is from the triple- 
gluon diagram. In contributions from the first two diagrams there are double-hard, soft-hard 
processes for each diagram. Here the second hard interaction refers to the quark-anti-quark 
annihilation into DY lepton pair and the initial interaction with the target gluon can be soft 
or hard, depending on the momentum fraction carried by the initial gluon. The first term 
has four contributions that include the interferences between the double-hard and soft-hard 
processes. For the second diagram, however, the soft-hard process is prohibited since the 
radiated gluon can only be induced by the first hard scattering. In this case, only double- 
hard process contributes. Notice that one can also obtain the double-hard contributions 
by considering two successive parton scatterings in a physical gauge as we demonstrate in 
Appendix O Using the above equation, we can get the leading twist-four contributions from 
annihilation-like processes to the DY differential cross section. 



z'^Tgg{x, xt) + Tjj^^ (x, xt) + (l - ^^i? i^^ ^t) , (26) 

where the nuclear twist-four matrix elements are defined as 

Tgq{x,xt) = J ^dy^dy^e'^^'+'^'^P^y' (l - e-'"*^^^^'"-^!")) (l - e-'"*^^^^ 

X l-{A\F^{y^)^PMl^My-)F^''{yT)\m-y2ny- - z/D, (27) 



T^¥ix,xt) = l^dy^dy^e'^P^y-e'^^^P^^y^-y^^ 



X ^{A\F^iy^)Moh^My')F^''iyi)\m-ymy- - vi), (28) 



-'^^^^X,Xt) = / ^rfyf d?/2"e^(^+^*)P''^ L-i^tP+y^ _2^-ixtP+{y -y,+y^) ^^-ixtp+{y -y,) 



X ^(A|F+(y2-)^,(0)7+^,(l/-)F+"(yr)|^)^(-2/2-)^(l/- " Z/D- (29) 

Here Tgq and Tg'q are related to the first and second diagram in Fig. HI respectively, and 

Tgq arises from the interference between them. Note that the interference between the soft- 
hard and double hard processes depend on the nuclear size and the momentum fraction Xt 
from either the target quark or gluon that is needed for the induced gluon radiation. The 
interference can be neglected when XiP'^Ra ^ 1 for large values of qt or short formation 
length for the gluon radiation relative to the nuclear size. In the small qt region or large 
formation time, however, the interference effect becomes important due to what is referred 



to as LPM effect |40]. The destructive interference is dictated by the gluon formation 



time, T/ = l/xtp^, which, in the coUinear limit, could become much larger than the relative 
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distance between the two scattering centers. Therefore the LPM effect will effectively cut off 
the transverse momentum of the virtual photon at ~ Q'^ /Ra 0] in the leading calculation 
of the DY lepton pair production. This will lead to an anomalous nuclear dependence of the 
DY cross section in h + A collisions. 

We also need to consider the interferences between single and triple scattering. Our cal- 
culations show that all the contributions from these processes are power suppressed because 
the leading term, i.e., the last term in Eq. f l22l) . vanishes. However, when combining all the 
contributions from both the double scattering and single-triple interferences, the first term 
in the coUinear expansion in Eq. f|T^ gives the eikonal contributions to the single scattering. 
Such contributions do not correspond to any physical scattering since they can be gauged 
away, but ensures the gauge invariance of the leading twist results. We list the full results 
of the single-triple interferences in the appendix. 

Similarly to the leading twist case, one also needs to consider the virtual contribution 
which corresponds to the quark self-energy correction. The results of virtual correction can 
be obtained by the requirement of unitarity. 
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(30) 



Including both radiative and virtual contributions, the final result for the annihilation-like 
process is given by 
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Here, the modified quark splitting function is 
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and 



AT - = 

91 



dz 

1 — z 
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(34) 



Notice that the infrared divergence is canceled between the real and virtual corrections in 
the same way as in the single scattering. 

The other twist-four processes that contribute to the modified beam quark distribution 
function are the Compton-like processes as shown in Fig. [51 The full results of these 
processes can be found in the appendix. Similarly to the calculation of annihilation-like 
processes, one can get the leading twist-4 contributions, 
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FIG. 5: Next-to-leading order Drell-Yan compton process at twist-four. 
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(35) 



which are analogous to the results from the annihilation-like processes. The first two terms 
correspond to the first two diagrams in Fig. [5l respectively. The third and fourth terms 
are the interference contributions between the two diagrams. The contribution from the last 
diagram in Fig. |5]is neglected as its contribution is power suppressed. Again, one can obtain 
the double-hard contributions by considering two successive single scatterings as shown in 
Appendix O Using the above results, for Compton-like process, Eq. (fTB|) finally reads 



dQ'^dx' Jo q^ '2tt Jx' ^ 

where the modified gluon splitting function is defined as 



-rfg/h{.i)Mg^qq{x' li, q^)Ho{x,p, q), (36) 
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with the twist-four matrix element 



'1 



z? + z' 



t: 



c_ 

99" 2 



2TcasCA 



(37) 



T. 



c 

91 



Z Tgq(^X, Xf) -\- Tgq iX-iXt) -\- zTgq (x,Xt). 



(38) 



Ca 

Diagrams in Fig. O contributing to each term in the above equation are similar to the 
annihilation-like processes. Here, the twist-four matrix elements are the same as in the 
annihilation-like processes as defined in Eq. (l27|) - (!29!) . In Compton-like processes, there is 
no virtual correction at 0{as), and the gluon splitting function is infrared finite. 



IV. MODIFIED QUARK DISTRIBUTION FUNCTION 

Summing up all the leading logarithmic contributions from single [Eq. (|H])] and double 
scattering [Eqs. (ET]) and fl5^ ]. we can define the medium-modified effective beam quark 
distribution function as 

2 A\ . / / 2n «s dq^ /-i d^ 



fq/h{x',fi ,A) = fq/h{x',fl ) + 7r — -r f<l/h{i)Mq^qg{x' li.qT) 

+ fgih{i)Alg^qq{x'/i,ql)\ , (39) 
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where fq/h{x',iJ, ) is the renormahzed twist-two beam quark distribution given in Eq. (|8]). 
The A-dependence of the medium-modified beam quark distribution function is imphcit 
through the medium-modified sphtting functions in the second term on the right-hand side. 
Assuming that multiple gluon bremsstrahlung induced by additional scattering in the nuclear 
medium can be resummed in the same way as in the vacuum, one can obtain the medium- 
modified DGLAP evolution equation for the beam quark distribution function 



dlrifi'^ 2n Jx' ^ 

where the modified splitting functions, 

%-^bc{z, /i^) = Ja~*bc{z) + A-fa-^bc{z, /i^), (41) 

are the sum of contributions from both the vacuum and the nuclear medium-induced 
bremsstrahlung. Notice that Eq. (1401) is formally the same as the DGLAP evolution equation 
in vacuum, Eq. ffTTl) . but its splitting functions have extra terms A7 which are caused by 
medium-induced gluon radiation. Such medium-induced terms will soften the beam quark 
distribution function as a consequences of radiative energy loss by the leading beam parton. 

As we mentioned before, the medium-modified splitting functions A7 depend on the twist- 
four matrix elements T4. In order to get a numerical estimate of the effect of initial parton 
energy loss, we need to know which is not calculable. Neglecting correlation between 
the quark and gluon distribution, one can assume the factorization of the twist-four matrix 
elements [3] in terms of quark and gluon distributions from two independent nucleons 
inside the nucleus, the later in turn can be related to quark transport parameter q in nuclear 
medium. Under these assumptions, the twist-four matrix elements can be expressed as 



Tgg{x,xt) = 2tx dy / d h dy^ pA^Vi ,b)pA{y ,b) 

J—oo Jo J—oo 

X sm\xtp^y^/2)f^/^{x, b) [xf^/^{x, b) + Xtf^/^ixt, b)] , (42) 



where, fa/N^^^^) parton distribution function per nucleon inside the nucleus. As 

illustrated in Fig. |6l 6 is the impact parameter of h + A collisions, yjf is the longitudinal 
position of the target nucleon where the soft gluon comes from and y~ is the longitudinal 
position of the target nucleon where the final quark-antiquark annihilation into DY lepton 
pair takes place. One has to integrate over the impact parameter b. The nucleon density 
distribution pA{y~,b) is normalized as 



00 



A= dy- d'bpAiy-,b). (43) 

J -00 Jo 

So far we have focused on the nuclear modification of beam quark distribution function 
in DY processes in h + A collisions due to multiple parton scattering and medium-induced 
parton energy loss. Multiple parton interaction and their coherence can also lead to nuclear 
modification of the target parton distributions. Gluon fusion in the evolution of the target 
parton distribution function can also lead to parton saturation [52|. Some of these effects 
are higher-twist as we have discussed so far and some are leading twist which are determined 
by the wave function of the nucleus and initial parton distribution functions. In this paper, 
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we will adopt a phenomenological approach for the nuclear modification of the target parton 
distributions and use the current parameterization of the nuclear parton distribution func- 
tions (nPDF) from global fit to the existing experimental data. In this case, the effective 
target parton distribution function per nucleon at impact parameter b is defined as 

/a%(^, b) = jfa/Aix, b) ^ R\{X, b)f,,^{x), (44) 

where fa/N{x) is the nucleon parton distribution function in vacuum and R\{x,b) is the 
nuclear modification factor for parton distribution functions inside a nucleus at impact 
parameter b. The impact parameter averaged parton distribution function in a nucleus is 
defined as 

fa/A{x) = fa/N{x) J d%t A{b) R'X{x , b) , (45) 

where t^(&) = / dy^ pA{y~ ,b) is the nuclear thickness function. Note that the nuclear mod- 
ification factor R'X{x, b) should contain the information about the isospin of the nucleus. 
We further define the generalized jet transport parameter as 



q{xt,y ,b) = ^ PA{y ,b)- xfgij^{x,b)\^^Q + XtfgiN{xt,b) 



« (46) 

assuming the limit Xf ^ 1, where go is the quark transport parameter in the center of 
the nucleus. With the above definition of quark transport parameter and nuclear quark 
distribution function, the twist-four matrix elements can be expressed as, 

6go/g-/jv(x) f ^ _ fV 

X sm\xtP+y^/2)R\ix,b). (47) 

In our numerical evaluation, we will use the Woods-Saxon distribution for the nucleon 
distribution 

PA{y,b)= P'^^ (48) 

1 + exp(^) 



Tgq{x,xt) ^ ^^ofqMx) I I dy j dy^pA{yi,b)pA{y ,b) 
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where r = a/5^~+~^, Ra ~ 1.1274^/^ fm is the radius of the nucleus, a ~ 0.537 the width of 
the nuclear skin and the constant po is fixed by normalization in Eq. f H5]) . 

Similarly, the other two matrix elements in Eqs. ( I28l) and ( l29l) can be approximated as 
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Pa (2/1 ,b)pAiy ,&) 



X sm' 



(xiP+yr/2)/?i(x,6). 



(50) 



Because of the LPM interference effect, the formation time for the induced parton radi- 
ation cannot be larger than the nuclear size. This will restrict the transverse momentum 
qt to a limited region ~ Q'^/A^^^. Below this limit, the contributions will be suppressed 
because of the destructive LPM interference. Such LPM effect in contributions that contain 
Tgq and Tg|'' give rise to an addition factor of A^^^ after integration over the transverse 
momentum. This nuclear enhancement leads to the quadratic nuclear size dependence of 
the parton energy loss. Relative to these two contributions and in the large nucleus limit, 
we can neglect contributions proportional to Tg^'^ which does not contain LPM effect [8|. 
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FIG. 7: Nuclear modification of the effective beam u-quark distribution functions 
fu{x',fi'^,A)/fu{x',ij,'^) versus x' in DY processes in p + yl collisions at Eiab = 800 GeV 
and scale /i = 4.5 (solid) and 2.5 GeV (dashed). The quark transport parameter is set at 
qo = 0.024 GeVVfmllJ. 



Shown in Fig. [7] are the effective nuclear modification factors of the beam u-quark distri- 
butions in p + W and p + Fe DY processes at Euh = 800 GeV. In this numerical calculation 
and in the rest of this paper, we will only consider the effective nuclear modified beam par- 
ton distributions due to double scattering as given in Eq. (15^ . The normal twist-2 quark 
distributions in nucleon are given by the CTEQ6L js^l parameterization, the nuclear modi- 
fication of parton distribution functions inside a nuclear target will be given by the EPS08 
43| parameterization without impact-parameter dependence. As shown in the figure, the 



nuclear modification of the beam quark distribution function is negligible in small x' region. 
The modification becomes significant in the large x' region for small factorization scale p, and 
large target nuclei. We also observe some small scale dependence of the nuclear modification. 
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Experimentally, the Drell-Yan differential cross section can be measured for different 
targets in hadron-nucleus collisions. The cross section ratios are analyzed for the study of 
nuclear effects, 



B(7^ BdapA^i+i- /dQ'^dx' 
Aa^ AdopB^i+i- / dQ'^dx' 



(51) 



In general, one would like to use the experimental data mp + p collisions as the reference in 
the denominator to study the nuclear effects inp + A collisions with a large nucleus target. 
One, however, often uses p + B with a light nucleus target B as an approximate of the p + p 
reference data (modulo isospin corrections) as in the Fermilab E866 [32] experiment. Under 
this approximation, the measured nuclear modification can be compared to 



B(T^ ^ Eg / dxfg/Ajx, iJ,^)fq/p{x', n"^, A)Ho{x, p, q) 
Aa^ AJ2q I dxfq/N{x, iJ,'^) fq/p{x' , iJ,'^)Ho{x,p, q) ' 



(52) 



where fq/p{x' , fi"^, A) is the modified beam quark distribution function defined in Eq. (j39ll . 
and fq/N{x,fi'^) is the normal twist-2 anti-quark distribution in nucleon (with the isospin of 
the nucleus B) which will be given by the CTEQ6L parameterization in our study here. 
As for the nuclear modification of parton distribution functions inside a nuclear target, we 
will use the EPS08 [43,] parameterization without impact-parameter dependence. 
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FIG. 8: Ratios of the DY cross section per nucleon in p + A collisions versus x' for Fe/Be 
and W/Be at Eiab = 800 GeV and invariant dilepton mass M = 4.5 GeV. The curves are 
theoretical results with (solid) and without (dashed) initial beam parton energy loss. The 
shaded bands correspond to % = 0.024 it 0.008 GeV^/fm [l5|. The experimental data are 
from the Fermilab experiment E866 32]. 



Shown in Fig. [8] are the calculated DY cross section ratios as a function of x' for different 
nuclei {Be, Fe and W) with (solid) and without (dashed) the effect of initial beam parton 
energy loss via multiple interaction, as compared to the Fermilab E866 experimental data 

. In our numerical evaluation of the effect of initial beam parton energy loss, the quark 



transport parameter go at the center of a large nucleus is the only free parameter and is fixed 
to go = 0.024 ± 0.008 GeV^/fm from the phenomenological study of parton energy loss and 
modified fragmentation function in nuclear DIS by the HERMES experiment 15|]. As shown 
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in the comparison between our calculation and the experimental data in the kinematic region 
of the E866 experiment, the dominant nuclear modification of the DY cross section is from 
nuclear shadowing of parton distribution functions inside large nuclei as given by the EPS08 



431] parameterization. For fixed invariant mass M, the fractional momentum x = M'^/x's 



carried by the target partons becomes smaller for large beam quark fractional momentum x', 
therefore strong nuclear shadowing of the quark distribution inside the target nucleus. The 
effect of medium-modified beam quark distribution caused by beam quark energy loss leads 
to further suppression of the DY cross section for large nuclei. However, with the quark 
transport parameter predetermined from the nuclear DIS experiment [l5|, the suppression 
due to initial beam quark energy loss is quantitatively small. This is consistent with other 
estimates of parton energy loss in DY process [33l^35|. The additional suppression only 



becomes considerable in large x' region in a large nucleus where one also see large nuclear 
modification of the effective beam quark distribution function as shown in Fig. [71 Since the 
parameterization of nPDF from global fitting included DY data, one should include the 
effect of beam parton energy loss in large x' or small x region. 

Without initial beam parton energy loss, the nuclear modification of the DY cross section 
as a function of the target parton momentum fraction x should more or less refiect that of the 
parton shadowing, which is shown in the left panel of Fig. [9l The slight x-scaling violation 
for different invariant DY dilepton mass M is caused mainly by the scale violation of the 
parton shadowing effect as given by the EPS08 [i^l parameterization. Parton shadowing 
becomes smaller for larger invariant mass due to the QCD evolution of the nuclear parton 
distributions which leads to less suppression of the DY cross section. When one takes into 
account of the beam parton energy loss, there is increased suppression of DY cross section 
at large x' or small x for fixed invariant mass M as shown in the right panel of Fig. [HI For 
fixed values of x, an increase of invariant mass M leads to an increase in x'. This should 
lead in turn an increased suppression due to beam parton energy loss according to Figs. [3 
Such increased suppression due to larger values of x' will counter the effect of scale violation 
of the shadowing (see the left panel). The final outcome is an approximate x-scaling of the 
suppression factor in small x region for different values of the DY dilepton mass. 

Nuclear shadowing of parton distributions in the target nucleus is significant only at small 
momentum fraction x. For fixed invariant mass of DY lepton pairs and moderately large 
beam parton momentum fraction x' or large target parton momentum fraction x, the effect 
of shadowing should be small. On the other hand the fractional parton energy loss will 



become larger for smaller beam parton energy x'Eiab [l5j. Therefore, at fixed DY dilepton 
mass M and moderately large beam parton momentum fraction x', the effect of the parton 
energy loss on DY cross section inp+A collisions will become more dominant for lower beam 
energy Eiab- This will enable one to disentangle the effect of initial-state energy loss effect 
from the nuclear shadowing. Shown in Fig. [10] are the predictions for the DY cross section 



ratios at Eiab = 120 GeV in the Fermilab's E906 experiment [5J, [55| with invariant dilepton 
mass M = 4.5 GeV. At this lower beam proton energy, the target parton momentum fraction 
X is large for moderately large beam parton fractional momentum x'. The effect of parton 
shadowing is indeed small as shown by the dashed line in Fig. [TUl On the other hand, the 
energy loss effect induced by multiple scattering is significant and the dominant cause for 
the DY suppression shown by the shaded bands in Fig. [TOl Therefore, the E906 experiment 
can provide an unambiguous measurement of the effect of initial state energy loss in DY 
cross section. This is consistent with another recent predictions of the energy loss effects at 
the E906 experiment j37| . 
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FIG. 9: Ratios of the DY cross section per nucleon in p+W and p+Be collisions versus x at 
Eiab = 800 GeV for invariant dilepton mass M = 4.5 (solid), 3.5 (dashed) and 2.5 (doted) 
GeV with only parton shadowing (left panel) and shadowing and initial beam parton 
energy loss (right panel). The quark transport parameter is set at % = 0.024 GeV^//m 
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FIG. 10: Predictions for the DY cross section ratios per nucleon in p + A collisions versus 
x' for Fe/Be and W/Be at {Eiab = 120 GeV) with invariant dilepton mass M = 4.5 GeV. 
The curves are theoretical results with (solid) and without (dashed) initial beam parton 
energy loss. The shaded bands correspond to % = 0.024 zt 0.008 GeV^/fm |15|. 



V. SUMMARY 

In this paper, we have calculated the differential cross section of the Drell-Yan dilepton 
production in hadron-nucleus collisions in the lowest order of pQCD, including the effect 
of initial multiple parton scattering and induced gluon radiation within the framework of 
generalized coUinear factorization. In nuclear medium, multiple parton scattering leads to 
energy loss of the beam parton which has a quadratic nuclear size dependence because of the 
LPM interference effect. Within the leading logarithmic approximation, we express the effect 
of parton energy loss in terms of an effective modified beam quark distribution function. The 
modification depends on the twist-four matrix elements which can be related to the normal 



19 



twist-2 parton distributions and jet transport parameter q. Using the values of q determined 



from nuclear DIS data by the HERMES experiment [15[, we found that the effect of initial 
beam parton energy loss is quite small for large invariant mass of DY lepton pairs m p + A 
collision at high-energies, except in the large x' region and for heavy nuclear targets. The 
Fermilab E866 experimental data are consistent with the effect of the nuclear shadowing of 
target parton distributions as given by the EPS08 parametrization. However, the effect of 
beam parton energy loss becomes significant for lower beam proton energy Eiab = 120 GeV 
at the Fermilab E906 experiment, especially at large x' and with heavy nuclear targets. For 
parameterization with improved accuracy of nuclear parton distribution functions through 
consistent global fitting of experiment data, one should also include the effect of initial beam 
parton energy loss in the DY process. 
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Appendix 

In this appendix we list the complete results for annihilation-like and Compton-like pro- 
cesses, we also provide an cross-check of the double-hard part through partonic quark-gluon 
and gluon-gluon scattering. 



Appendix A: complete results in annihilation-like process 



There are total 9 diagrams for central-cut, and 7 diagrams for left- and right-cut each. 
The central-cut diagrams are shown in Fig. [Hand their contributions are: 
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All the left-cut diagrams in annihilation- like process are shown in Fig. [TT] and their 
contributions are: 
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where, 



H.,, = -s-erne,-^^-^^^!.. (A26) 

7^^5 = e'("^+"*)P^^"e'("^-"*)P''^i"e-'"^P''^^"0(-y2")^(y2" -l/D; (A27) 
- 2 2Ca 1 + ^' ql + kr-qT 8(27r)% 

i^AL6 = a.aeme — -— — A28 

2 1-2; g2 (gy + fcy)2 

(A30) 

7^^7 = _e^(^B+xt)p+s/-g-i(xf+x0p+yre^^^P+y2"^(_y-)^(l/2" -l/i")- (A31) 



Tf 2 2^ 1 + 1 8(27r)% 

-nAL7 — Oigaem.eqt'A- 2 AT2 — 



= -^{2qT ■ kr - k^), = -qi ' ^ ' 

The right-cut diagrams are the complex conjugate of left-cut diagrams shown in Fig. [TTl 
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their contributions are: 

J^ARl — Olg'^em.egL'F- 2 7^2 — 



X 



AR2 



alttf^mel ( Ci? ^ 



Ca\ 1 + gf^ - (1 - 2)A;t ■ Qt 8(2 



s-en.., V^^ 2 J 1-z ql[qr - (1 - 
1ar2 = e'^^^+^'^P*y-e^''^^^y~^'y'^^e{y- -y^)e{y:[ -y:,)- 



'■R2, 



HaR3 

Tars = e 



X _ g*(a;t-a;'i3-a;£;)p+{j/i -j/j ) 

2 2<^A 1 + gr + -zfcT • Qt 8(27r)=^ - 

tt^ttemeo ^ o/^ r \o Ar2 



Jxtp+{y2 -y ) 



X e 

Har4 = culaemelCF 



I _ ^i(xt-xn+XF)p'^iyi -1/2 ) 



1 + ^2 1 8(27r) 



'RA, 



2 ^2 r_ _ - Q^o^^2 



H 



AR5 



2 2<^a 1 + ^' QT-kT-qT 8(27r) 



2 1 - Z (g^ - iVc' 



H 



AR6 



2 a'^A 1 + QT + kT-qT 8(27r)^ - 

r>i r,,^ r*<„ ato ^ RG) 



2 1-z g2(^^ + fc^)2 iV2 

T7 2 2^ 1 + ^' 1 8(27r)% 

HaRI — Ois'^em^q'^A- 2 M2~^'^^ 

i — Z Qrp iV^ 

J^^7 = -e^(^B+xt)p+y-^-^xpp+iy--y^)^.xtp+iy^-y-)0f^y~ _y-^0l^y~ _y-- 



(A33) 

(A34) 
(A35) 

(A36) 



(A37) 
(A38) 

(A39) 
(A40) 
(A41) 
(A42) 

(A43) 
(A44) 

(A45) 
(A46) 

(A47) 



Appendix B: complete results in Compton-like process 

There are same amount of diagrams in Compton-hke and annihilation-hke processes, a 
total of 9 diagrams for central-cut, and 7 diagrams for left- and right-cut each. The central- 
cut diagrams are shown in Fig. |5]and their contributions are: 



Hcii 


= a^a 


Hc22 


= a^a 


Hc33 


= ala 



Ca 


2z^- 


2z + l 


B(27r)2 




1 {qr- 


zkrY 




Cf 


2^2 _ 


2z + l 


B(27r)2 




1 {qx- 






Cf 


2^2 _ 


2z + l 


B(27r)2 



qx 



Nr. 



Ic22, 



(Bl) 
(B2) 
(B3) 
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Hci2 


2 2 


1 




-2z + 1) 


Hci3 


_ 2 2^ A 


1 




-2z + l) 


Hc21 


_ 2 2^ A 


1 

c 


-(2.^ 


-2^ + 1) 



[g|-(l + ^)fcT-gT + ^4]8(27r)^ ^^^^^ (B4) 
ql - zkr ■ qr 8(27r)% 



- (1 + z)kT ■ Qt + -2^T 8(27r)^y 



/ci3, (B5) 
(B6) 



—' ^ - C21 

{qT — kxY^qT — zkxY 
-rj 2 2(n '^A ^ m 2 o , ,^ qT-kT-qT 8(27r)% 

^C3i = «^«..eJ^-^(2.^ - 2. + l) tr'%ty ^^3i, (B8) 



and 



(BIO) 



The left-cut diagrams in Compton-like processes are shown in Fig. [121 Their contribu- 
tions are: 

HcLi = ala,mel^^^{2z' -2z + 1)^^^!/^,, (Bll) 



0{-y2)0{y2-yiy, (B12) 

(B13) 



L2; 



Tr 2 2 '^A /o 2 o , n - (1 - ^)^T " QT 8(27r)2 

"(7f2 = o^zafrnez—, — 2^: — 2^: + 1 1 

' ''2(iV2-l)^ ^g2[g-^_(l_^)fc^]2 

X ^-ixtv*y2 ^i{^E-xt+XD)p^ {Vi -y2) _ ]^ ■ (B14) 

T7 2 2 /o 2 o , 1 N 4 + ■ gr 8(27r)2 

^^-3 = ",".™e,^p^(2. - 2z + (B16) 

Ici3 = e'<«+""'''*»"e-«'»'*'"."-»J)e(-j,-)e(y--y-) 



3i(-ZF-Xt+a::D)p+(j/^ -y2 ) _ I 



(B16) 



HcLA = ala,^el—^{2z^ -2z + 1 - A^/^^, B17 

IcLi = -e'("«+"*)^'^^"e'("^-"*)P^^i"e-*"^P^^2-^(_^-)^(^- _^-); (b18) 

2 2£f-CV2.. 2 o , ,, qT-kT-qT 8(27r)% 

i^CLs = a.aeme,— —(22; -22; + l)— (B19) 

~ [qr - kr) qT 

7c^5 = e'(-s+-')p''s^"e*(-^-')P+2^re-*"^^'^^2-^(_^-)^(^- _^-); (B20) 

7f -^^r. g2CWV2 o ql + kr-qT m^Yj 

^CL6 - a^aemeg — r^^ — [2z - 2z + L)—— — i^g, lB21j 

- J- qHqT + krr 
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FIG. 12: Left-cut diagrams in compton-like process. 



'■CL6 



2 



HcLi = a^aemeq—^ -{2z -2z + l) 



1 8(27r)2 
— ^7 — -'^7, 



'CLl 



(B22) 
(B23) 
(B24) 



The right-cut diagrams are the complex conjugate of left-cut diagrams shown in Fig. 
and the results are: 



T7 2 2 

ncRi = a^aemeq 

ICRI = 



-{2z^ -2z + l) 



1 8(27r) 



X 



1-e 



ixtp+{y^ -J/2 )p«2;tp+(y2 -y ) 



-^Ci?2 



-(2^2 



2 _ 2^ + l) gT-(l-^)fcT-9T8(27r)^ - 



^ ^ixtp+(y2 -y ) 



X _ gi(^t-XB-XE)p+iyi -1/2 ) 



(B25) 

(B26) 
(B27) 

(B28) 



(B29) 
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Hem 
Icm 

Hem 
Icm 

Hem 
Icm 

Hem 

Icm 
Hem 

Icm 



C 



A 



"s«-e,2(iV2-l) 



^iixB+xt)p+y-^-ixf^P+{y--y-)Q(^y- _ y-)e{y- - y^) 
\ _ QKxt-XD+XF)p+ {y^ -y^) 



m, 



y. ^ixtp+iy^ -y ) 



2 2Cf-C^a/2,^^2 , in #-^r-gT 8(27rl 



-{2z' -2z + l) 



'm, 



^i{xB+xt)p+y-^ixEP+(y^-y^)^ixtp+{y^-y-)Q^y- _ y-)e(jj- - y^); 



q^ + kr- qr 8(27r)% 



Im, 



gi(xs+xOp+y-g-ix^p+{y--y-)gix,p+(s,--y-)^^^- _ y-^0(^y- _ ^-). 



2 2 



Cf 



'^m - 1 



(2^2 - 2^ + 1; 



1 8(2 



qr 



,i(xs+a'0p+2/ g-i^FP+K -J/2 )ei^tP+(y2 )0(^y- _y-^0(^y- _y^y^ 



(B30) 

(B31) 

(B32) 

(B33) 

(B34) 

(B35) 

(B36) 

(B37) 
(B38) 

(B39) 



Appendix C: Cross-check of double hard processes 

One can extract the double-hard subprocesses in annihilation-hke and Compton-hke pro- 
cess from Eq. ( l26l) and Eq. ( l36l) . 



dQ'^dx' 



dQ'^dx' 



dq^ as d^ 



X 



J dxHo{x,p,q) 
\Caz + Cf(1 - zf Thh{x, xt 



q^ 271 J X' ^ 



1 + ^2 2nas 



Z Nr. 



/ dxHo{x,p,q) 



dq^ as di 



gj. 2tx Jx' ^ 



'-fg/hiO 



2'Kas 

m-1 



X, Xf) 



(CI) 



(C2) 



where. 



THH(a;,a;,) = J ^dy^dy^e'^'^'-y- e'^^^^^^y^-y^^ 

X ^(A|F+(y2-)^,(0)7-^^,(y-)F+"(yr)|A)^^(-y2")^^(y- - yr 



(C3) 



/ dy^pA{yN)fq/A{,x)xtfg/N{,Xt) 



with paIvn) being the nucleon density distribution. A factorized form of gluon-quark cor- 
relation is assumed in the above equation. With this assumption, one can factorized the 
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double-hard subprocess as, 



J dxd^Ho{x,p,q)fg/h{Ofg/Aix) J dy^pAiVN] 



1 -\- z 

X na1xtfg/N{xt)- 

i — z 



2 r 



dql 
dz— 

Qt 



da 



S(0) 

hA^l+l- 



dQ^^ 



dQ^ 



dyj^PA{yN)daqN, 

J2 I dxd^Ho{x,p,q)fg/hiOfg/A{x) J dy^pAiVN) 

Cp Ca 



(C4) 



X iralxtfg/Nixt) + (1 - zf 
dyJjpA{yN)d(ygN. 



m - 1 m-i 



z{l-z) 



dz 



dq^ 



"''^hA^l+l- 



dQ"' 



(C5) 



Here da^A^j* / dQ"^ is the Born cross section defined in Eq. ([6]). daqN and dagN represent 
the cross section of the parton-nucleon scattering. Since the gluon is physical and has finite 
momentum fraction Xtp'^ in the double-hard process, the higher-twist results have simple 
and intuitive partonic interpretation as stated in the above equations. This comes from 
the assumption that the twist-4 contributions from double scattering can be factorized into 
two successive and isolated single scatterings. In the following, we rederive the twist-4 
contributions from double-hard scattering based on the above assumption 49 . 
Consider the parton-nucleon scattering 



a(p') + N{p) ^ c(g) + dik) + X 



where, 



p' =[0,p'-,0], p=[p+,0,0], q 
The parton nucleon cross section can be written as 



2zp'- 



zp ,qT 



daaN = I daagfg/N{x)dx 

f 
2s 



dxfg/N{x) — \Mlg^^a{s,t, 



u 



(27r)32go 



2tt6[{p' + xp - qY 



9 



dz 



j^Xtfg/N{xt)\M\^g_^^^{s,t,u)^ ^^ dq'j,, 



(C6) 
(C7) 

(C8) 
(C9) 

(CIO) 



where Xt = qp/l'^P ' ~ ^)] because of the on-shell condition of the final parton. The 
Mandelstam variables in this partonic process can be recast as 



s = {p + xp)' 



Qt 



— t = {q-p'Y = - — , u = {q-xpf = — . 

z[l — z) z 1 — z 



(Cll) 
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In the annihilation-like and Compton-like double scattering processes, the first hard scat- 
tering is the quark-gluon Compton scattering and gluon-gluon fusion, respectively. The 
partonic matrix element for these two processes are 



|M| 



s, t, u) 



\M\l^^^{s,t,u) 



u 



Nc us 

CaI + z^ CfI + z^ 
+ 



N, (1 - zf N, z 



Cf t' + u' 



Ca f + u" 



m - 1 tu 



Cf z^ + U-z? 



C, 



(C12) 



zl-z) 



m - 1 



\z' + (1 - ^) 



(C13) 



Substituting the above equations into Eq. (ICSp . one can obtain the quark-nucleon and 
gluon-nucleon cross section. 



Ca , Cf ,^ 



1 + z"^ , dq^ 



da, 



gN 



]xtfg/N{xt) 



C, 



Ca 



1 - z 



-z{l - z) 



-dz- 



4 ' 

(It 



z^ + zf 



(C14) 
(C15) 



which are equivalent to the results in Eq. (]C4p and Eq. (ICSp . 
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